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The scaled distribution of the smallest eigenvalue in the Laguerre orthogonal and symplec- 
tic ensembles is evaluated in terms of a Painleve V transcendent. This same Painleve V 
transcendent is known from the work of Tracy and Widom, where it has been shown to 
specify the scaled distribution of the smallest eigenvalue in the Laguerre unitary ensemble. 
The starting point for our calculation is the scaled fc-point distribution of every odd labelled 
eigenvalue in two superimposed Laguerre orthogonal ensembles. 

1 Introduction 

A quite old result in mathematical statistics concerns the eigenvalue distribution of random 
matrices of the form A = X T X where X is a n x N (n > N) rectangular matrix with real 



entries. First it was proved by Wishart [27] that with 

n N 

(dX) := J] II dx i^ (dA) := J[ da jk 

3=1 k=l l<j<k<N 

denoting the product of differentials of the independent elements, the change of variables from 
the elements of X to the elements of A takes place according to the formula 

/ \ (n-JV-l)/2 

(dX) = (detAJ (dA). (1.1) 

From this a description in terms of eigenvalues and eigenvectors can be undertaken by introduc- 
ing the spectral decomposition 

A = OAO T 

where the columns of O consist of the normalized eigenvectors of A, and A is the diagonal matrix 
of eigenvalues. About a decade after the work of Wishart, it was proved by a number of authors 
(see e.g. 0) that 

N 

(dA)= J] lAfc-Aj-l JJ dA i(° T d°)- U- 2 ) 

l<j<k<N j=l 



A significant qualitative feature of (1.2) is that the eigenvalue dependence factors from that of 
the eigenvectors. 

Suppose now the elements of X are identical, independently distributed standard Gaussian 
random variables so that the corresponding probability measure is proportional to 



n N N 

J] II e~^ k/ \dX) = e~^ xTx \dX) = e~^^(dX). (1.3) 

j=l k=l 
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Noting that det A = Yh=i > substituting (|1.2| ) in ( |L~l| ) and then substituting the resulting 
formula for (dX) in gives the now standard |l7| result that the eigenvalue probability 

density function (p.d.f.) of the matrix A = X T X is given by 



1 

C 



N 

n*. 



(n-JV-l)/2 Xj/2 



n \^ 

l<j<k<N 



A 



3\i 



(1.4) 



where C denotes the normalization and Xj > (j = 1, . . . ,N). This is referred to as the real 
Wishart distribution, or alternatively as the Laguerre orthogonal ensemble (LOEn), the latter 
name originating from the occurence of the classical Laguerre weight function X a e~ x (up to 
scaling of A) and the invariance of ( |1.4| ) under the mapping A i— ► OAO T . 

Another random matrix structure which leads to the p.d.f. (|1.4j) is the block matrix 




(1.5) 



It is straightforward to verify that in general (1.5) has n—N zero eigenvalues, while the remaining 
2N eigenvalues come in ± pairs. It is similarly easy to verify that with X distributed according 
to ( |1.3| ) the positive eigenvalues are distributed according to ( |1.4[ ) but with Aj i— > X? and an 
additional factor of HjLi \Xj\. Hence the square roots of the positive eigenvalues of ( |1.5| ) are 
distributed according to (ITT 



Over the past decade the p.d.f. ( |1.4| ) has found application in a number of physical problems. 
One example is in the theory of quantum transport through disordered wires, in which the 
matrix product X T X for X a N x N random matrix modelling the top right hand block of the 
transmission matrix occurs in the Landauer formula for the conductance (see e.g. Another 
is in quantum chromodynamics, where the structure (|l.5|) models a random Dirac operator in 



the chiral guage, and the number of zero eigenvalues is prescribed [26|. 

Because A is positive definite and so has positive eigenvalues the eigenvalues near Xj = 
are said to be near the hard edge. For eigenvalues in this neighbourhood, it is known that with 
n — N fixed, the statistical properties tend to well defined limits in the N — > oo scaled limits, 
where the scaling is 

Xj h+ Aj /AN (1.6) 
I?], ||]. In fact the scaled fc-point distribution function is known exactly [12]. Thus with 



K h (X,Y) 

D\{x,y) 
I\{x,y) 

fi(x,y) 



(n- N - l)/2 

J a (A 1 /2)yi/2j/(yi/2 ) _ x 1 /2j^(A 1 /2)j a (yi/2) 



2(X - Y) 



(1.7) 



p\° ) Eh (x 1 ,...,x k ;a) 



J S^x, z) dz - -sgn(x - y) 

S\(x,y) I\{x,y) 
D\{x,y) S\(y,x) 

v ( 1 \ k LOEjv/ Xl X k 



we have 



P[k) Eh ( x i' ■■■ ,x k ](a- l)/2) = qdet [fi(xj,x k )]j,i=i,...,k 



(1.8) 
(1.9) 
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where qdet is the quaternion determinant introduced into random matrix theory by Dyson 
and the superscipt "h" denotes the hard edge scaling (1.6). 



In this work we will compute the exact distribution of the smallest eigenvalue in the scaled 
LOE as specified by the Appoint distribution Ql.S| ) in terms of a certain Painleve V transcendent. 
We will also compute the same distribution for the scaled Laguerre symplectic ensemble (LSE), 
which before scaling is specified by the eigenvalue p.d.f. 

j=l l<j<k<N 

This p.d.f. results from positive definite matrices A = X*X when the matrix X has the 2x2 
block structure 

(z w \ 
—w z J 

of a real quaternion (each eigenvalue is then doubly degenerate as well as occuring in db pairs). 
The explict form of the corresponding scaled /c-point distribution function is given in ]l^]; in 
taking the scaled limit with scaling ( |1.6| ) the ensemble LSEjv/2 is considered rather than LSE^v- 

Crucial to our study is knowledge, in terms of a Painleve V transcendent, of the distribution 
of the smallest eigenvalue in the scaled Laguerre unitary ensemble (LUE). Before scaling, the 
latter distribution is specified by the eigenvalue p.d.f. 

1 N 

-nA?e- Aj n i**-^i 8 . (i-io) 

J=l l<j<k<N 

and results from positive definite matrices A = X'X when the matrix X has complex elements. 
The corresponding scaled /c-point distribution function 

P(k) (xi,...,x fe )- j bm ) ^— ) P[k) ^— j 

has the explicit form || 

pg Eh (xi, . . . , x k ) = det[K h ( Xj , Xi)]j,l=i,...fc. (1-11) 



where K h is given by (1.7). 



In general the probability that there are no eigenvalues in an interval J, E(0; J), can be 
written in terms of the corresponding /c-point distribution by 

00 (-l) k f f 
E(0;J) = 1 + ^^—^ / da*... / dx kP(k) (x u ...,x k ). (1.12) 

k=i k - Jj Jj 

For p(u\ a k x k determinant with entries g(xj,x k ) the structure ( |l,12j ) is just the expansion of 
the Fredholm integral operator on J with kernel g(x,y). Thus it follows from ( |1.11| ) that 

E%(0; (0,s);o) = det(l — K h ) 

where E^O; (0, s); a) denotes the probability there are no eigenvalues in the interval (0, s) for the 
scaled LUE (the subscript 2 characterizes the LUE via the exponent on the product of differences 
in ( 1.1C| )), while K h denotes the integral operator on (0, s) with kernel K h (x, y). With p™ m (s; a) 
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denoting the distribution of the smallest eigenvalue in the scaled LOE (/? = 1), LUE {(3 = 2) or 
LSE {(3 = 4) we have in general 



d -*h, 



p^(s;a) = --E$(0;(0,s);a) 

so to compute p y g m {s; a) it suffices to compute E^(0; (0, s); a). 

For general a > — 1, ^(fy (O) s )i°0 has been computed in terms of a Painleve transcendent 



by Tracy and Widom [22]. The Painleve transcendent is denoted by q^ (in [£2|] the subscript h 
is not present), and specified as the solution of the nonlinear equation 

s(q 2 h - !)(«£)' = q h (sq' h ) 2 + ±(s - a 2 )q h + ±*gg(gg - 2) (1.13) 



subject to the boundary condition 



Qh(s) ~ 7^ rs a/2 . 

yhV ',_»o+ 2 a r 1 + a) 



That is a Painleve transcendent follows from the transformation [22] 

/ \ i + y(x) 2 

Qh{s) = - rT , s = x 

l-y(x) 

from which one can deduce that y(x) satisfies the Painleve V equation 



= (7r + i )(y)--y+ 2 — {ay + -)+ — + 

\2y \ — y' x x z \ x y — 1 

with a = — /? = a 2 /8, 7 = and (5 = —2. The result of [22] is that the Painleve transcendent % 
specifies £3 via the formula 

E^(0;(0,s);a) = e W (-^ J\\og S /t)(q h (t)) 2 dt). (1.14) 

In this work we will show that E\ can also be evaluated in terms of qh(t)- Specifically, we 
obtain the formula 

^(0; (0, s); (a - l)/2)f = £ 2 h (0; (0, s); a) exp ( - \ jT ^ dt) . (1.15) 

With E\ and known in terms of the probability £4(0; (0, s);a) can also be expressed 
in terms of qh(t) by using the formula |l3| 



E{(0; (0, .); a + 1) = 1(^(0; (0, .);<„- l,/2) + -^g&g^— ) . (1. 16 ) 

Thus 

(El(0; (0, S ); a + 1)) 2 = £ 2 h (0; (0, a); a) cosh 2 (± jT ^ dt) . (1.17) 

Here -E^ 1 is computed by scaling the ensemble LSE7V/2 according to ( |l.6|) . Also, with E\(l\ (0, s); a) 
denoting the probability that there is exactly one eigenvalue in the interval (0, s) for the scaled 
LOE, we have the inter-relationship [JDJ 

E\{1; (0, s); (a - l)/2) = £ 4 h (0; (0, a); a + 1) - ^(0; (0, a); (a - l)/2). 
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Substituting ( 1.15 ) and ( 1.17 ) shows 

(X h (l; (0, s); (a - l)/2)) 2 = £ 2 h (0; (0, s); a) sinh 2 (± jf* dt) . (1.18) 

Crucial to our derivation of (|1.15|) is a reworking of the derivation of Tracy and Widom 
[23] giving the probability Ef(0; (s, oo)) in terms of a Painleve transcendent. Here Ef(0; (s, oo)) 
denotes the probability that there are no eigenvalues in the interval (s, oo) for the scaled Gaussian 
orthogonal ensemble (GOE). As the density falls off rapidly as s increases, the region (s,oo) is 
said to be a soft edge, thus explaining the use of the superscript "s" in Ef. 

The ensemble GOEjy refers to the eigenvalue p.d.f. 



i N 

3=1 ±<j<k<N 



This is realized by N x N real symmetric matrices, with diagonal elements having the Gaussian 
distribution N[0, 1], and independent off diagonal elements having the distribution N[0, l/y/2]. 
Tracy and Widom begin with the quaternion determinant expression [16] for the Appoint dis- 



tribution function in the finite GOE (we also draw attention to the recent work |2[ in which 
Painleve type recurrence equations are obtained for the analogue of the probabilities Ef and E\ 
in the finite system). Instead, inspired by the observation of Baik and Rains || that the square 
of the distribution of the largest eigenvalue in GOE is equal to the distribution of the largest 
eigenvalue in two independent, appropriately superimposed GOE's, technically the ensemble 

even(GOE 7V UGOE 7V ), (1.20) 



we take as our starting point the /c-point distribution of ( 1.20 ), scaled at the spectrum edge. 

Now the /c-point distribution of ( |1.20| ) is an ordinary determinant, whereas the /c-point dis- 
tribution of the GOE is a quaternion determinant. Furthermore the elements of the determinant 
contains terms familiar from the analysis of E% given in pi]] ; this is also true of the quaternion 
determinant but the former involves only a subset of the latter. These facts together provide a 
simplified evaluation of E\. The power of this derivation is demonstrated by its application to 
the evaluation of E\. We find that each step used in the derivation of E\ has an analogous step 
in the case of E\ and this leads to ( |1.15 ). 



We begin in Section 2 by providing the evaluation of the /c-point distribution for the ensemble 
( 1.20| ) in the scaled limit at the right hand soft edge, as well as that for the ensemble 



odd(LOEjv U LOE^v) (1.21) 

in the scaled limit at the hard edge. In Section 3 we begin by using the evaluation of pnA 
obtained in Section 2 as the starting point for the evaluation of E\(Q\ (s, oo)), and then proceed 
to mimick this calculation to evaluate £^(0; (0, s); a). In Section 4 our evaluations ( 1.1 5| ) and 



( 1.17 ) are related to previously known results. 



2 The ensembles even/odd(OEiv(/) U OEjy(/)) for / classical 

Let OEjv(/) denote the matrix ensemble with eigenvalue p.d.f. 

N 

n/te) n ixk-xji (2.i) 
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We see from (|1.4j ) and ( 1.19 ) that the LOE is of this form with 

/( x ) = x a e~ x/2 , (x>0, a := (n-N- l)/2) (2.2) 
while the GOE is of this form with 

fix) = e~* 2 l\ 

In fact the four special choices of / 



e 



x / 2 , Gaussian 



... , , x^ a x )/ 2 e x / 2 (a: > 0), Laguerre , . 

■ /( - r) i (l-x)( a - 1 )/ 2 (l + 2 ;)( b - 1 )/ 2 (-l<x<l), Jacobi [ } 

(1 + x 2 )-( Q+1 )/ 2 , Cauchy 

(note that the exponent of x in the Laguerre case has been renormalized relative to fl2.2|)) have 



been shown in 13] to possess special properties in regards to the superimposed ensembles 

even(OE^(/)uOE iV (/)) and odd(OEjv(/) U OE N (f)) (2.4) 

(amongst other superimposed ensembles). In particular the /c-point distribution is given by a 
determinant formula with the same general structure in each case. 

To present the formula for the fc-point distributions, some additional theory from 13] must 
be recalled. In particular, it is found that each of the weight functions (|2.3j ) is one member of 
a pair (/, g) which naturally occur in the study of the superimposed ensembles. Explicitly the 
weight functions g are 

e~ x2 , Gaussian 



x a e x [x > 0), Laguerre 
(1 - x) a (l + x) b (-1 < x < 1), Jacobi 
(l + x 2 )- a , Cauchy, 



Now, let {p n (x)}n=o,i,... denote the set of monic orthogonal polynomials associated with a par- 
ticular weight function g, and let (p n ,Pn)2 denote the corresponding normalization. Then it is 
shown in that for the ensemble even(OE7v(/) U OEjv(/)) 



■N-2 



. T-r , , , . pi{xj)pi(xj) pjy_i(gj).Fjv_i(gj) 



(2.6) 

i,j=l,...,n 



where 



F JV _ 1 (x)= f; (pllzk^x), I_{ x ):= f M [ X f(t)dt. (2.7) 



Similarly for the ensemble odd(OEjv(/) U OEjv(/)) the /c-point distribution is again given by 
the formula fl2,6|) but with J_ in (|2.7|) replaced by 



I + ( X ):=M r f(t)dt. (2.8) 
9{x) Jx 

The summation in ( |2.6D can be evaluated according to the Christoffel-Darboux formula, and 
the corresponding scaling limits are well known ||. To compute the scaled limit of the quantity 
-F/v_i in (2.6), we first make use of results from the work [Q] to provide the explicit evaluation 
of the coefficients 

/■oo rx 

(p,,J_) a := / dxf(x)pi(x) / dtf(t), 
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applicable in all the classical cases (p.q). First we note 



f(t)dt 



1 



sgn(x - t)f(t) dt + - f(t) dt =: 4> (x) + s , 



which allows us to write 

(Pl,l-h = 
Now results in Q give 

4>o (x) 

where 



f(x)pi(x)4>o(x) dx + s / f(x)pi(x)dx. 



7k(Pk,Pkh = ' 



1 

2 ' 



Gaussian 
Laguerre 



i(2fc + 2 + a + 6), Jacobi 



a — 1 — k, Cauchy 
This allows us to immediately evaluate the first term in ( p-10| ). 

It remains to compute the second term in (2.10). With the notation 



81 := 2 



f(x)pi(x)dx. 



this term is given by 2s"oS;. Consider first the case I even. With 

1 



(j)l(x) :-- 



sgn(x-t)pi(t)f{t)dt 



we know from flll that 



4>2k(x) 



9{x) 



En 



72fc ri/ fe =i(72«-i/72«) f( x ) t^hl=i ^ (P2i/+i,PaH-i)a 



72Z-1\ P2v+l{x) 



(2.9) 



(2.10) 



(2.11) 



Forming the ratio (j)2k{x) I §2k-2(x) an d taking the limit x —* oo shows that S2fc/*2fc-2 
l2k-2llik-\ an d thus we have the evaluation 



% = *o n — 



1 sO) P2fc(a;) 



To evaluate si, I odd, we recall the formula JTJ 

72fc 72k J{X) {P2k,P2k)2 

Taking the limit x — > oo implies s*2fc+i = (72^-1 /72fc)*2fc-i an d since s_i := this gives 

%+l = 0. 

Thus the second term in (2.10) is fully determined. 
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Substituting the evaluation of (pi, I-)2 obtained from the above working in (2/7) shows 



771 / \ V"^ IT|=1 (72/-l/72j) / \ I o ~2 \- n j= o(72j/72j+l) , , , . 

F N -i(x) = — ^ — P2u+\(x) + 2s ^ — =£-, p2i(x). (2.12) 

^° v=[(N-l)/2] \P2v+l,P2v+l)2 l=[N/2] \P2UV2l)2 

From this result we read off that 

, f ^n!=r 2)/2 (TSU-i/TSa) ^even 
^-l,fjv-l 2= 7 ° 9 w n/2 (2-13) 

{ 2 S inl=r 1)/2 (72,--2/72i-i) ivodd 

so all quantities in the expression ( |2.6|) for are now known explicitly. 

In the case of the ensemble odd(OE7v(/) U OEjy(/)) the definition (2/7) of Fjv-i has J_ 
replaced by 1+. Noting 

) 1 ^00 1 rco 

f(t) dt = -- / sgn(x - t)f(t) dt+- f(t) dt = -Mx) + s 

£ J — CO ^ J — CO 

which differs from ( |2.9| ) only in the sign of the first term, we see by revising the working which 
led from ( |2.9| ) to ( |2.12j) the only modification needed to the formula ( 2.12| ) is that a minus sign 
be placed in front of the first term (and similarly in ( |2.13| )). 

2.1 Guassian ensemble at the soft edge 

In the Gaussian case 

f(x) = e-* 2 ' 2 , g{x)=e- x \ Pl (x) = T l H^x) , (p,, Pi ) 2 = tt^HI, (2.14) 
where H[(x) denotes the Hermite polynomial. The soft edge scaling is || 

X 



x = ( 2iV ) V2 + ^§1/6 (2 - 15) 



so we seek to compute 



where E denotes the ensemble (|1.2C| ) and the r.h.s. is given by (|2.6| ) with the substitutions ( p. 14 ). 



Regarding the summation in (^1]) , we know from the study of the GUE at the soft edge that 

i 

Oo^Nm (p IjPI )2 - (2N)1 £ +x '% N l% - K[X > Y) > 

1=0 \rnri)<!. y =(2N) l l 2 +Y/2 l l 2 N 1 / b 

where, with Ai(x) denoting the Airy function 

:= Ai(x)Ai-(,)-Ai-(,)Ai(,)^ 
x-y 

This is obtained using the Christoffel-Darboux summation formula and the asymptotic expansion 

e- x2 ' 2 H n {x) = 7r- 3 / 4 2"/ 2 + 1 / 4 (n!) 1 / 2 n- 1 / 12 (vrAi(-n) + 0(n- 2 / 3 )) (2.17) 
where x = (2n) 1 / 2 - ii/2 1 / 2 n 1 / 6 . 
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It remains to compute the scaled limit of the term involving Fpf—i in (2.6). Now substituting 
the values of 7^ and (pk,Pk)2 from ( 2.1 1] ) and Q2.14 ), and noting Sq = ir/2, ( 2 . 1 2| ) reads 



00 u\ 00 1 

F N -l(x) = 7T 1/2 (2 U + ^y H 2u+l{x) + 7T ^ ——^(a;). 



^=[(AT-l)/2] 



(2.18) 



i=[AT/2] 



Next we want to combine this result with ( 2.13| ). For definiteness take N to be even. Then we 
see that 



(g(x)g(y)) 



/2 PN _ 1 {x)F N _ l {y) 
(PN-ii F N _i) 2 



e- x2 l 2 2^ N ^H N ^{x)(A 1 {y) + A 2 



(2.19) 



with 



A 2 (y) :-- 



-2/2/2 00 



(N/2 - 1)1 ^ (N - 1 + 2i/)! 

1 



t/=0 
l/2 p -2/ 2 /2 00 



7T ' e 



(JV/2 - 1)! ^ 2^ v + 2 '(A r /2 + /)! 

We remark that the summation defining Ai(y) (with the lower terminal u = replaced by v = 1) 
occurs in the study of the soft edge distribution at f3 = 1 [12], and furthermore a procedure has 
been given to compute its asymptotic behaviour with the scaling ( 2.15 ), the key ingredient of 
which is the asymptotic expansion ( |2.17| ). 

For the x-dependent terms in ( |2.19[ ), ( |2.17| ) gives 



(v), 
H N +2i(y)- 



2.20) 



e - 2 /2 2 -(N-D HN _ l{x) 



~ 7T 



a: =(2Ar)i/2 + x/2 1 /2Ari/6 
For Ai(y), first note that for large N 

(N/2-l + u)l(N -l)! 1 / 2 



l/4 2 -(AT-l)/2+l/4 (iv _ 1)! l/2 iV -l/12 Ai(x) _ 



(N - l + 2u)\ 1 / 2 (N/2- 1)! 



Then use of ( p. 17 ) with n = N + 2u — 1, —u ~ Y — 2v/N 1 ^ shows that the sum becomes the 
Riemann approximation to a definite integral, and thus 



(N - l)! 1/2 ^i((2iV) 1 / 2 + Y/2 1 / 2 N 1 ' & ) ~ 2^' 
Hence 



-l)/2 2 l/4 -1/4^-1/12^^ 

2 



Ai(l" - v) dv. 



e -^/2 2 - ( N-l) HN _ i{x)My) 



x=(2N) 1 / 2 +X/2 1 / 2 N 1 / 6 r> 
a={ 2 JV )l/2 + y /2 l/2 Jv l/6 



-Ai(X) [°° Ai(Y - v) dv. 
- Jo 

(2.21) 



Similarly, for A 2 (y), noting that for large iV 

(iV + 21)! 1 / 2 (iV-1)! 1 / 2 1 _ t 

2 N + 2l (N/2 + l)\ (N/2 - 1)1 ~ (27)V2 2 ' 

and using ( |2.17|) with n = N + 2u, — u ~ Y - 2v/N 1 l' i we find 

(iV - l)! 1 / 2 yl 2 ((2Ar) 1 / 2 + W2 1 / 2 iV 1 /6) ^ J_ Al u 2 N)^ 2 + Y /2 l ' 2 N 1 ^) 
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and hence 



lim 



1 



tf_oo 2V2AT1/6 



e -^/2 2 -(N-l) HN _ i{x)My) 



1 f 00 

a; =(2iv)i/2 +x/2 i/2 Jv i/6 = -Ai(X) / Ai(y-t;)du. 

y=(2iV) 1 /2 + y / 2l/2 J vl/6 Z •'0 



The contributions ( |2.2l| ) and ( |2.22| ) thus reinforce, so after adding to (|2.16| ) we obtain 



p^ ) OE)2s (X 1 ,...,X k ) = det 



(K^X^Xt) + Ai {Xj) J™ Ai(A ; 



v) dv 



j,l=l,...,k 



(2.22) 



(2.23) 



2.2 Laguerre ensemble at hard edge 

In the Laguerre case 

f(x) = x {a - 1)/2 e- x/2 , g(x) = x a e- x , (x > 0) 
p,(x) = (-l) l l\Lf (x), (p l ,p l ) 2 =T(l + l)T(a + l + l), (2.24) 
where Lf denotes the Laguerre polynomial. At the hard edge the appropriate scaling is specified 



Xu 



(2.25) 



by (|1.6|) , so the task is to compute 

(LOE) 2 h,^ v , ,. / 1 \ k E f X l 

where E denotes the ensemble ( 1.21 ) and pi^\ on the r.h.s. is specified by Q2.6| ). 

The scaled limit of the summation in ( p.6| ) at the hard edge occurs in the study of the 
LUE and is known 
asymptotic formula 



Thus making use of the Christoffel-Darboux formula and the large n 

e -y/2 y a/2 L a {y) „ ^2 J^^y (2 . 26 ) 



one finds 



N-2 



lim — — 



Pi(x)pi(y) 



N-*ooAN f^ Q (pi,pi) 2 



x=X/4N 
y=Y/4N 



K h (X,Y) 



(2.27) 



where K h is specified by (|Q|). 

The first step in computing the term involving fV— 1 in (|2.6| ) is to substitute the formulas 
( 2.24| ) i n (|2.12[ ) modified so that there is a minus sign before the first term (recall the paragraph 
below (p.iap . This gives 

„ . . 2a! ~ 2 2 »{a/2 + v)\v\ Ta 

F N _i{x) = j-j^ -r-M^ — -rr^-LSL-u (x) 



|2Q ((a-l)/2)! 2 (a/2)! « 

l=[2V/2] 



(2/)! 



a l 



2 2l l\(a/2 + l)\ 



We see from this formula and (2.13) (for definiteness in the latter we take N to be even; recall 
that in this case a minus sign must be inserted) that 

(JV-l)! 



v / Cpjv-i,-Pjv-i)2 



2 N - 2 ((N - 2)/2)!(o/2 + (JV - 2)/2)\ 



2 2v (a/2 + v)\v\ 



E 



(a + 2^ + 1)! 



(5(y)) 1/2 ^+i(y) 



+2 



i/=(JV-2)/2 

^((a-l)^)! 2 ^)! 2 * 



a! 



1 2 



E 



(20! 



l=N/2 



2 2l l\{a/2 + l)\ 



(g(y)) 1/2 L a 2l (y) 



{g{x))y 2 L a N _ 1 {x) 



(N-l)\ 



2 N - 2 ((N - 2)/2)!(o/2 + (iV - 2)/2) 



■(Bi(y)+B2(y)), (2.28) 
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where B\ and B2 denote the two terms in the final brackets of the line before. 
Consider the function B\{y). From Stirling's formula we have 

2 2v {a/2 + v)\v\ 



(a + 2v + l)\ 



7T 



l/2 2 - (0+1)^-0/2- 1/2 



Using this and ( |2.26| ) we that the sum defining B\(y) is the Riemann approximation to a definite 
integral and we find 



£1 



Y 



TT 



1/2 /N\l/2 



7T 1 / 2 /iV\l/2 1 f 00 



v4iV/ 2 a / 2+1 V 2 
Another application of Stirling's formula shows 

(N-l)\ 



2 a/2 V 2 / yV2 



/■oo 

/ J a (i)dt. (2.29) 



2 (a+l)/2 7r -l/2 Ar (l-a)/2 



2^- 2 ((iV - 2)/2)!(a/2 + (N - 2)/2)! 
so we see from further use of fl2T26l ) together with (|2.29|) that 

nm ^^^(5( x )) 1 ^ 2 47V^7V-l( X )2 A '- :J ((iV-2)/2)!(a/2+(Af-2)/2)! ^M?/) 

A similar analysis applied to B2(y) shows 

lim A r-+oo ^(5(^)) 1/2 ^7V-l(^) 2 A'-2(( A r_2)/2)!(a/2+(Af-2)/2)! ^(2/) 



x=X/AN 
y=Y/4N 



(2.30) 



a: = X/4JV 
i/ = y/4JV 



which thus reinforces ( 2.30| ). Thus adding twice ( 2.30 ) to ( 2.27[) we have 



pf k ^ E)2h (X 1 ,...,X k ;(a-l)/2) 



det 



(k\x j ,x 1 ) + Ja{ ^ 



2^1 J^xi 



Ja(t)dt 



j,l=l,...,k 



(2.31) 



(2.32) 



3 Gap probabilities at the spectrum edge 
3.1 The probability El (0; (s, 00)) 

It was remarked in the Introduction that the probability E\ (0; (s, 00)) has been computed in 
terms of a Painleve II transcendent by Tracy and Widom P3|| . Expliclity let g s denote the 
solution of the particular Painleve II differential equation 



(q s )" = sq s + 2(q s 



(3.1) 



subject to the boundary condition q s (s) ~ Ai(s) as s — > 00. Then it is shown in [ 23 1 that 

(£f(0; (s, 00))) 2 = £!(0; (s, 00)) exp ( - J™ q s (t) dt) (3.2) 



where [21 



El(0;(s,oo)) =exp(- / (t-s)q 2 s (t) 

J S 



dt 



(3.3) 
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Here we will use the evaluation of the /c-point distribution ( |2.23| ) to provide a simplified derivation 
of (|3.2|) while still following the essential strategy of ]23| . 



By definition of (1.20) it follows that 

, 2 



( E GOE N{Q . ( a>oo)) y =i *ven(CSOE w UGOB w ) (0;(S)Oo))j 
which with the scaling ( p,15| ) implies 

(i5?(0; (s, oo))) 2 = i^( GOE ) 2s (0; (s, oo)). (3.4) 

Now, recalling the determinant formula ( p. 23 ), we see from ( [1,12 ) and the text immediately 
below that £( GOE ) s can be written as the determinant of a Fredholm integral operator. Thus 

(^?(0; (a, oo))) 2 = det (l - (if 8 + A 5)) (3.5) 



where X s is the integral operator on (s, oo) with kernel fl2.16|) while A is the operator which 
multiplies by Ai(x), while B is the integral operator with kernel Ai(y — v) dv. 

Removing (1 — K s ) as a factor from (|3.5|) and recalling [tO] 



det(l -if 3 ) = £f(0;(s,oo)) 

we obtain 

(^(0;(s,oo))) 2 = E|(0;(s,cx>))det (l - (1 -if 8 )- 1 ^® B) 

= E s 2 (0;(s,oo)){l- J^(l- K s )- 1 A[y]B(y)dy) (3.6) 

where the second equality follows from the fact that (1 — if s ) _1 A[y] is the eigenfunction of the 
operator (1 — K S )~ 1 A tg) B, so the eigenvalue is 

(l-^-^BCi/)^ 

Analogous to the notation of |2l| we put 

<A s (x) := A(x) = Ai(x), Q s (x) := (1 - if s ) _1 7l[x] 

so that 

/•oo /*oo ry 

/ (l-^ s )- 1 A[y] J B(y)dy= / dyQ s (y) cf> s (v) dv =: u s e (3.7) 



(the notation u 8 - without the superscript s - is used for an analogous quantity in p3| ). Note 
from fl3.6|) that with the notation fl3.7|) we have 



(^(0; (a, oo))' = ^(0; (0; (s, oo))(l - <)■ (3.8) 

Following 1 23], our objective is to derive coupled differential equations for u e and the quantity 

ry 

dyp s (s,y) / f(v)dv (3.9) 

J— oo 

where p s (x,y) is the kernel of the operator (1 — X s )" 1 . These equations will involve 

/oo 
dyp s {s,yW(y), (3.10) 
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which in is shown to be the Painleve II transcendent specified by the solution of (3.1), and 
their derivation relies on the formula [21] 



9 Q S (V) = -qs(6 + (y -s) + p s (s, y)) , (3.11) 



ds 

where 5 + {y — s) is such that 

/oo 
S + (y-s)f(y)dy = f(s), 

as well as the formula [Bill 



I- + I- + Plfffav) = -Q s (x)Q s (y). (3.12) 



- ds dx dy > 
Now, differentiating ([T?]) with respect to s we have 

' <P\v) dv + J dy(-Q s (y)) J ^ <f>\v) dv = -q s q s e (3.13) 



where to obtain the second equality use has been make of ( |3.11| ) and the definition ([?!]). We 
now seek a formula for (ql)'- Making use of ( |3 . 1 2| ) in ( |3.9| ) shows 

/•oo fi ry rco ry 

(<£)' = - dy—p s (s,y) ct> s (v)dv-q s dyQ s (y) ct> s (v)dv 



OC 



dy 

p s (s,yW(y)dy - q s u s t 
= q s (l-<) (3.14) 
where the final equality follows from the definitions (3/7) and ( |3.10|) . 



As q s is known, the system of equations ( |3.13| ) and ( |3.14| ) fully determines u s e and q s e once 
boundary conditions are specified. Now Q s (y) is smooth, so we see from ( |3.7| ) that 



u s e ~ (3.15) 

as s — > oo. On the other hand p(s,y) = 5 + (s — y) + R(s,y) where R(s,y) is smooth, so for 
s — ► oo 



OC 



OO 



/oo 
Ai(«)du = l. (3.16) 
-oo 



The unique solution of the coupled equations ( 3.13; ) and ( |3. 14 ) satisfying ( |3.15| ) and ( |3.16| ) is 
easily shown to be 

u\ = l-e-^% q s e = e~^ (3.17) 

where 

poo 

fjL s := / q s {x)dx. 

J s 

Substituting the evaluation of u\ from ( |3.17 ) in ( |3.8|) reclaims (|3.2j), as desired. 



3.2 The probability E%(0; (0, s); a) 

All the steps leading to the rederivation of (|3.2| ) given in the previous section have analogues for 
the probability (0; (0, s); a) which lead to the evaluation ( 1.15| ), 

First, the analogue of ( |3.4[) is 



E*(0; (0, S );a)) 2 = ii;( LOE ) 2h (0;(0, S );a), 
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while use of the determinant formula (|2.32| ) in ( 1.12 ) then gives 



(^(O; (0, s); (a - l)/2)) 2 = det (l - (K h + C ® D 



(3.18) 



Here iT h is the integral operator on (0, s) with kernel ( |1.7|) , while C is the operator which 
multiplies by J a (y / y), while D is the integral operator with kernel 



Ja(t)dt. 



(3.19) 



Recalling (2 



det(l - K h ) = E%(0; (0,a);o) 
we see that analogous to (|3.6|), (|3.18| ) can be rewritten 

(E?(0; (0, s); (a - l)/2)) 2 = £ 2 h (0; (0, a);a)(l-J\l- K h y l C[y]D(y) dy) . (3.20) 

Analogous to the notation of |2^] we put 

h (^) := C(x) = J a (v^), Q h (2/) := (1 - K h )~ l C[y]. 

After changing variables t = ^Ju in ( |3.19| ) we see that in terms of this notation 

1 ( s 1 f°° 1 

(1 - K )~ 1 C[y]D(y) dy = - dy Q h (y)— / du ^\u) =: t& (3.21) 

4 Jo \jy Jy y/U 

and in turn this latter notation used in ( 3.20[ ) gives 

(tfftO; (0, s);(a- l)/2)) 2 = £*(0; (0, s); a)(l - (3.22) 

Q h (s) =:g h = f dyp h ( S ,y)0 h (y), 
Jo 

which in |22j is shown to be the Painleve V transcendent specified by the nonlinear equation 
( 1.13j) , analogous to ( 3. 11| ) we have \22\ 

§-Q h (y) = g^ + (y - s) + p h (s,y)). 

Use of this formula in ( |3.21 ) then shows 

1 f s 1 f°° 1 

(*£)' = T«h<£ q^-= dyp h (s,y)— du^{u) (3.23) 

which is the analogue of ( |3.13| ). 

Next we seek a formula for the derivative with respect to s of q^. For this purpose we note 
from [^] that 

x^- P h (x,y) + s^-p\x,y) = -±( yp \ x ,y)) + Iq\ x )Q\ v ) 



Now, with 



dx 



ds 



dy 
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(c.f. (|3.12| )). This formula applied to ( |3.23D shows 

d 



s(g, 



r s / d \ l f 00 l 

/ (— (^(s,?/) — / du^4> h (u)+q h v^ 
Jo ^dy ) yjy J y yju 

1 rs 1 roo 1 rs 

- / dyp h (s,y)— / dn^0 h (n)- / dy p h (s, y)0 h (y) + g h u 

2 Jo V2/ Jy V « JO 



l q * - qh (l - u h € ) 



(3.24) 



The coupled equations ( 3.23|) and ( p. 24 ) must be solved subject to the s — > boundary conditions 

roo 1 yoo 

^ ~ 0, ~ / -n=4>{u) du = 2 J a (v) dv = 2. (3.25) 

jo V u Jo 



The occurence of yfsq^ in ( 3.25 ) suggests we introduce 

£ ■■= ^ 

in ( |3.23|) and ( 0-24 ). Doing this gives the system of equations 



Introducing the new independent variable 



s 1 

Yj^ q h (x) dx 



o x 



we see that (3.26) reduces to the system with constant coefficients 



d h 1 ~h ^ ~h / -i h\ 



The solution satisfying ( |3.25| ) is 



(3.26) 



(3.27) 



(3.28) 



The stated formula (|L15D for E^(0; (0, s); (a - l)/2) now follows by substituting this evaluation 
of in (|3~22|). 



4 Discussion 

4.1 Special values of a 

Edelman @ was the first person to obtain the exact evaluation of E\ (0; (0, s); a), albeit for two 
special values of a only, namely a = —1/2 and a = 0. In terms of the scaling ( |1.6[ ) the results of 
§ are 

£i h (0;(0,s);-±) = e~^ s+ ^ (4.1) 
^(0;(0,a);0) = e" s / 8 . (4.2) 
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Subsequently it was shown by the present author (|| that E^(0; (0, s); a) for a G Z>o can be 
expressed as a 2a-dimensional integral. Explicitly 



E 



2a 

it 

' 3=1 



s 1 / 2 cos 



l<j<£:<2a 



2a 
(4.3) 



where 



2a 



r(3/2)r(3/2+j72) 

r(i + i/2) 



We remark that well known integration procedures (see e.g. [^4|]) allow this integral to be ex- 
pressed as a Pfaffian. Such Pfaffian formulas, deduced in a different way, have been given in 

hi. 



The formula ([LID relates E^(0; (0,s); (a - l)/2) to £#(0; (0,s);a), so it is appropriate to 
consider the evaluation of the latter for special values of o. Analogous to ( |4.3| ) we have that for 
a G Z> H 



S 2 h (0;(0, S );a) = e- s / 4 



fJ_) a i_ / TT e ^ 2 cos^ TT le^-e^^...^ (4.4) 



This integral can easily be written as a Toeplitz determinant, with the result [11] 

£#(0; (0, s); a) = e~ s/i det k_ fe (s 1/2 ) 



j,fc=l,...,a 



(4.5) 



where I n (x) denotes the Bessel function of purely imaginary argument. As an aside, it is 
interesting to note that the integral (^4|) is the generating function for the enumeration of 
various combinatorial objects, including quantities related to random permutations [19], random 
words [25] and random walks [10]. 

The formula ( |1.14| ) gives 



4(4 lo ^<°;<M;")) - ( 



(4.6) 



so ( [4.5[) implies that for a G 2>o, q\ can be expressed in terms of Bessel functions. The simplest 
case is a = 0, when we have 



£ 2 h (0;(0,s);a; 



-s/4 



(4.7) 



Substituting in (|4.6|) gives []22| 

9h(s) = 1 (4.8) 
and substituting this in ( |1.15| ) we reclaim fl4.ip . The next simplest case is a = 1 when we have 

E%(0;(0, S );l)=e- s / 4 Io( S 1/2 ) (4.9) 

and so 



1-4- 



d / d 



ds V (is 



^log/o^ 1 / 2 )). 



For this to be consistent with ( 1.15 ) we see that the identity 
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must hold. This in turn is equivalent to the statement that y := log Iq{s 1 I 2 ) satisfies the nonlinear 
equation 

Asy" + As(y'f + Ay' - 1 = 0, 
a fact which is readily verified using Bessel function identities. 

Special evaluations are also known for E^(0; (0, s);a) in the case a 6 Z>o 0. These evalua- 
tions are in terms of a certain generalized hyper geometric function based on Jack polynomials, 
while in the case a even Pfaffian formulas are also known [18]. In quoting from these results, 
one must be aware that E\ is defined starting with the ensemble an d scaling according 

to flOP . This means that the results of M require some rescaling of s. Doing this, we note from 
H that the simplest cases are a = and a = 1, when we have 



£4(0; (0,s);0) 
£i h (0; (0,s);l) 



-s/8 



Mb 



s 
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e s / 8 cosh ■ 



(4.10) 
(4.11) 



The result ( [4.10 ) can only be related to ( 1.17| ) in the limit a — > — 1~, since for a < — 1 
£?2 (0; (0, s); a) = 0. However, as the limiting forms of the quantities on the r.h.s. are not known 
we cannot readily check the consistency with ( [4.1C ) . In contrast the consistency between ( |1.17 ) 



and ( 4.11 ) is immediate upon recalling (|4.7| ) and ( 



ED 



4.2 Connection between Ep and E^ 



In previous articles |^, 12] it has been noted that for a — ► 00, after appropriate rescaling of the 
coordinates, the scaled /c-point distribution function for the infinite Laguerre ensemble at the 
hard edge becomes equal to the scaled /c-point distribution function for the infinite Gaussian 
ensemble at the soft edge. Note that in the symplectic case the scalings are done starting with 
the ensembles GSE^r/2 and LSEjv/ 2 - Let 

(a-l)/2, = 1 
a(P) = { a, (3 = 2 

a + 1, j3 = A 

Explicitly, it was checked that the scaled /c-point distribution function for the infinite Laguerre 
ensemble at the hard edge, with parameter a 1— ► a((3) and after the rescaling of coordinates 

x^a 2 - 2a(a/2) 1/3 x, 

equals the soft edge distribution functions for the corresponding Gaussian ensemble results in 
the a — ► 00 limit. We must therefore have 



lim E%(0; (0, a z - 2a(a/2) 1/3 s); a(3) ) = E 8 B (0; (0, 00)). 



(4.12) 



To verify ( 4.12| ), we first recall some additional results from [22]. Write 



(W) 



A 1/2 



so that, after integrating by parts, ( |1.14 ) reads 

£§(0;(0,a);o) =exp(-^ R h (t)dt 



(4.13) 



(4.14) 
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Then a := sRh(s) is shown to satisfy the particular Painleve III equation in a form (for an 
account of the latter see fllSl |5| ) 



2/ '\2 



Similarly let 



so that (3J3) reads 



(sa") 2 + <x'(<7 - sa'){Aa' - 1) - a\a 



,(*) = (-i?^)) 1/2 



£1(0; (s,oo)) = exp ( - i? s (t) 



0. 



It is shown in [21] that i? s satisfies the particular Painleve II equation in a form 

(<) 2 + 4W(i^ 2 -^ + i*s) =0. 



(4.15) 
(4.16) 
(4.17) 

(4.18) 



Substituting ( f4.14j ) and ( [4 . 1 7| ) in ( ^.12| ) with /3 = 2 we deduce that the validity of the latter is 
equivalent to the statement 



2a(a/2) 1/3 R h (a 2 - 2a(a/2) i/3 t) ~ RJt). 

a — *oc 

This can be verified by introducing the function 

? ( S ) = 2a ^l )V \ (a 2 - 2a{a/2) 1 '' i s) ~ 2a(a/2) 1 / 3 R h (a 2 - 2a(a/2) 1 / 3 s 



1/3, 



(4.19) 



into ( 4.15[ ) and taking the limit a — > oo. One finds the differential equation ( [4.18 ) results with 
Rs = a{s). 

In terms of ( f4,13| ), the evaluations ( |1.15| ) and ( |1.17 ) read 



(^(0;(0 )S );(a-l)/2)) 2 = £ 2 h (0; (0, s); a) exp ( 



((tR h (t)) 



nl/2 



(^(0;(0,s);a + l)) 2 = ^(0; (0, s); a) cosh 2 ( - / """ 'J. ' r// ). 



\/t 

1 /- (ti? h ^))') 1/2 
2io " 



Making use of ( 4.12f ) in the case (3 = 2 (which has just been verified), and ( |4.19| ) together with 
( 4, 16] ), we see that 



Inn ^(0;(0,a 2 - 2a(a/2) 1 / 3 s); (a - 1)/2)J = £§(0; (s, oo)) exp ( - y ft(i)cft 
a lirn (^(0;(0,a 2 - 2a(a/2) 1/3 s); a + l)) 2 = E|(0; (s, oo)) cosh 2 (- J°° q B (t) tit). 
We recognize the right hand sides in these expressions as E%(0; (s, oo)) for (3 = 1 and 4 respec- 



tively ^] (recall eq. (|3.2|) ; as noted in [ 13 1 £?|(0; (s, oo)) can be deduced from Ef (0; (s, oo)) and 
#1(0; (s, oo)) because of the validity of an inter-relationship analogous to ( |1.16| )). 
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